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ABSTRACT 

In Aristotelian Electrodynamics (AE), due to radiation over- damping, the velocity rather than 
the acceleration of a charge is determined by the local electromagnetic field. Treating electrons 
and positrons separately, AE seems to give a faithful description of the flow of charges in a pulsar 
magnetosphere. AE might allow to calculate the pulsar emission in full detail, at all energies 
except the radio. 

But AE is not a full theory - the plasma production processes must be added. Here we show 
that the pulsar magnetosphere and the resulting emission are sensitive to the plasma production 
rate near the light cylinder. In strong pulsars (high production rate, due to photon-photon 
collisions), the number of field lines closing beyond the light cylinder decreases. This reduces 
both the pulsar efficiency (the ratio of the pulsed bolomctric luminosity to the spin-down power) 
and the characteristic photon energy, in overall agreement with the Fermi data. 

For weak pulsars (low production rate), our numerical implementation of AE works stably 
and seems to be ready to calculate the lightcurves and spectra. But for strong pulsars, the code 
(included) needs improvement or better understanding. 



1. Introduction 

To "solve the pulsar", one needs to describe 
the production, motion, and radiation of electrons 
and positrons. Aristotelian Electrodynamics (AE) 
describes the motion and curvature radiation of 
charges in a strong electromagnetic field. AE is 
applicable to pulsars. 

AE appears to be ready (after going 3D, high 
resolution) to calculate the lightcurves and spec- 
tra of weak pulsars (insignificant pair production 
near the light cylinder). For strong pulsars (high 
rates of pair production near the light cylinder), 
one needs to add the relevant QED kinetics. This 
seems to be doable (§4.3), but there is another 
problem. At high production rates, our AE code 
behaves in an unclear manner. Since the code was 
written by a numerical amateur, there must be 
room for improvement. 

In §2 we describe Aristotelian Electrodynamics, 
in §3 we use AE to calculate the pulsar magneto- 
sphere, in §4 we discuss how AE can be used to 



calculate the pulsar emission. 

2. Aristotelian Electrodynamics 

In pulsar-strength electromagnetic field, be- 
cause of strong radiation damping, velocities 
rather than accelerations of positrons and elec- 
trons are given by the local electromagnetic field: 

E x B ± {B B + E E) 

v± = WTWl • (1) 

Here the scalar E and the pseudoscalar B are 
the proper electric and magnetic fields defined by 

Bl - El =B 2 - E 2 , B Q E Q = B • E, E > 0. (2) 

Equation (fTJ) gives unit velocities (speed of light 
c = 1). The charges actually move slightly slower 
than light (see §4). 



'-More precisely: (i) where Eq is large, the charges move 
at high Lorentz factors, (ii) where Eo is small, the sign of 
Bo frequently changes and the charges move at an average 
velocity below c. 
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Equation (JT|) says that in the frames where E 
is parallel to B, the charges move at the speed 
of light parallel and antiparallel to E. As they 
move along ((TJ, the charges are pulled by E 
and dragged by the curvature radiation. We 
will assume that a charge reaches the terminal 
Lorentz factor before it moves over the character- 
istic length scale of the field. For pulsars, where 
all the action is near the light cylinder (see below) , 
this requirement of instantaneous parallel dynam- 
ics reads 



, R 

L sd > L e ( — 



2/3 



(3) 



where L s d is the spin-down power, R = £ is the 
light cylinder radius, ft is the angular velocity of 
the star, r e — -S-j- = 2.8 x 10 _13 cm is the classical 
electron radius, L e = — 8.7 x 10 16 erg/s is the 
"classical electron luminosity". In astrophysical 
notation, the applicability condition @ reads 



L 34 » 5.7 x 1Q-°R 



5p2/3 



(4) 



where L34 is the spin-down power in units of 
10 34 erg/s and P ms is the pulsar period in ms. The 
applicability condition (fj| is satisfied, and by a 
large margin, by all pulsars in the Fermi catalog 
(Abdo et al 2010). 

3. Using AE to describe pulsars 
3.1. AE equations 

To calculate the pulsar magnetosphere, one 
needs to solve Maxwell equations. Maxwell equa- 
tions can be solved (numerically) if one knows the 
electric current j. According to AE, the current is 



J = - p-v-, 



(5) 



where p + is the positron charge density, /?_ is the 
absolute value of the electron charge density, and 
v± is given by ([T]). 

The charge densities should be calculated from 
the continuity equations 

P± + v-(p±v±) = g, (6) 

where Q is the plasma production rate. 

3.2. Previous AE calculation of pulsars 

In our previous AE calculation of the pulsar 
(Gruzinov 2012a), instead of modeling the plasma 



production Q, we postulated an expression for the 
charge- normalized number density po, 



Po 



P- 



(7) 



(We used p = yj ' p 2 + f 2 (B 2 + E 2 )/r 2 , where 
p = p + — p- is the electric charge density, r is 
the spherical radius, and / is a dimensionless pa- 
rameter which we varied.) 

Once the plasma density p is postulated, equa- 
tion ^ becomes an Ohm's law - an expression 
giving j in terms of the electromagnetic field - 



pE x B + Po (B B + E -E) 



B 2 



El 



p = V • E. 



(8) 



(9) 



Now the magnetosphere can be calculated in terms 
of the electromagnetic degrees of freedom only. 

The result of the calculation was non-trivial. 
It turned out, that the magnetosphere hardly de- 
pends on the assumed "fiducial multiplicity" / . 
In particular, we claimed to have calculated the 
universal value of the pulsar efficiency (the ratio 
of the pulsed bolometric luminosity to the spin- 
down power): 



0.5 



1 + 5 sin 2 6 ' 



(10) 



where 9 is the spin-dipole angle. The result- 
ing 0-averaged efficiency e ~ 0.1 is in an over- 
all agreement with the data, especially if one ig- 
nores brighter pulsars. Since in AE the drag comes 
from the curvature radiation, one can also give 
a crude estimate of the photon energy (see §4) 

E - cihe-iL] d R-i - SL^P^lGeV. This again 
roughly agrees with the data if the brighter pulsars 
are ignored. 



3.3. Why the postulated multiplicity model 
is wrong 

Given the insensitivity to the assumed multi- 
plicity, one might think that the actual plasma 
production rate is irrelevant. But it is easy to see 
that this cannot be right. 

Suppose the plasma production rate near the 
light cylinder Q is high, 



Q > flpcj, 



(11) 
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where pcj ~ B/R is the Goldreich- Julian (1969) 
density, B is the magnetic field at the light cylin- 
der. What is the fate of all these constantly added 
charges? 

The charges cannot be removed faster than 
light, and the plasma density po will be much 
higher than pcj- The proper electric field Eq has 
to be small in most of the volume - otherwise the 
damping rate q = j • E = cpqEq would be much 
greater than L s d/R 3 - At small Eq, and consid- 
ering the axisymmctric case for simplicity, we sec 
from (fTJ) that the poloidal motion of the charges 
is along the poloidal magnetic field. Then at least 
one of the species has to move into the equatorial 
current layer (see Fig.l). Once inside the current 
layer, this species will have to experience the full 
potential drop of the current layer. This again 
results in more damping power than L sc i- 

The above argument is not in contradiction 
with the model of (Gruzinov 2012b), where we 
show that certain plasma production rates Q do 
give a fixed-multiplicity AE. These models require 
pair annihilation - negative Q values must be al- 
lowed. Then both species move into the equatorial 
current layer and get annihilated there. But in real 
pulsars the annihilation rate should be very small, 
and such models should not be relevant. 
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Fig. 1. — Weak pulsar - no plasma production 
near the light cylinder. The star rotates at half 
the speed of light. Thin black - poloidal current 
isolines. Thicker yellow - poloidal magnetic field. 
Thick red - electric potential. 



3.4. New AE calculation of pulsars 

For the aligned rotator, we have numerically 
integrated the full axisymmetric AE system: (fTJ 
03 H]) plus Maxwell equations 



B = -V x E, 

E = V x B - j. 
The plasma production was 

Q = ocxQx + a 2 Q2, 



(12) 
(13) 

(14) 



where Q\ models the avalanche plasma production 
near the star (Ruderman and Sutherland 1975) 
and Q2 models the plasma production by photon- 
photon collisions near the light cylinder. 
For no particular reason we used 



Qi = 0(1.4r fl 



r)6{E -E c )(l 
2 ~ 1 + (ilrf ' 



(15) 
(16) 




Fig. 2. — Stronger pulsar - significant plasma pro- 
duction near the light cylinder. 
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where r is the spherical radius, r s is the radius 
of the star, E c is the critical field for the near- 
star avalanche, is the angular velocity of the 
star. Further details of the model are given in the 
Appendix. 

With no plasma production near the light cylin- 
der (Fig.l), we get a magnetosphere similar to 
(Gruzinov 2012b). Since the equatorial current 
experiences a large potential drop, this pulsar will 
have a high efficiency. 

For significant plasma production near the light 
cylinder (Fig. 2), the number of field lines closing 
beyond the light cylinder decreases. The equato- 
rial current experiences a smaller potential drop, 
this pulsar will have a lower efficiency. The cal- 
culated damping rate, J d 3 r (p + + p-)E , indeed 
drops by about a factor of 2. 

At low plasma production, the simulation re- 
sults seem to converge just as nicely as in a 
purely electromagnetic simulation of (Gruzinov 
2012ab). But at significant plasma production, 
our code probably needs improvement or better 
understanding. The convergence of the code (with 
time step and resolution) is unclear to the author. 
The results of the calculation given in Fig. 2 have 
the status of an illustration. But since the un- 
derlying physics seems to be right, we decided to 
show this illustration. 

How small can the pulsar efficiency become at 
high multiplicity? A crude estimate of §4 suggests 
that efficiency e can't decrease too much below the 
ideal value of e ~ 0.1. Observationally, even the 
strongest pulsar (Crab) has efficiency e <~ 0.01 - 
just a factor of 10 below the mean ideal value (X- 
ray dominated; the light cylinder of Crab can be 
opaque to high-energy gamma-rays). 

From a purely theoretical perspective, the min- 
imal efficiency problem is unclear to the au- 
thor. Our results seem to indicate that at high 
plasma production rates the AE magnetosphere 
approaches the standard zero-efficiency force-free 
magnetosphere of Contopoulos et al (1999) and 
Spitkovsky (2006). But one needs to remember 
that the "periodic axisymmetric pulsar" (Gruzi- 
nov 2011) has a 100% efficiency even at an arbi- 
trarily high multiplicity. 



4. Pulsar emission 

According to AE, the Poynting energy is ex- 
tracted at the current layer just beyond the light 
cylinder. AE predicts the terminal Lorentz factor 
of the charges 7: 

ceE ° = l^ 4 ' ( 17 ) 

where R c is the radius of curvature of the tra- 
jectory of a charge. R c is also predicted by AE. 
Knowing 7 and R c , one calculates the spectrum of 
the primary gamma-rays (usually <~ GeV ). 

In weak pulsars, the primary gamma-rays must 
be dominating the emission, while being irrelevant 
for the plasma production - the plasma is pro- 
duced only by the near-star avalanche. In strong 
pulsars, the situation seems to be more compli- 
cated. 

4.1. Weak pulsars 

The Fermi catalog lists 6 pulsars with P ms <~ 
300 and i 34 ~ 1. These have the photon cutoff 
energy ~Gev. For 3 of these pulsars, the data of 
Possenti et al (2002) gives a crude upper bound 
for the X-ray luminosity. Neglecting the geomet- 
ric factors, one then gets a crude estimate for the 
plasma production near the light cylinder. Com- 
paring to the Goldrcich- Julian density, we con- 
clude that these pulsars are weak. 

These pulsars still have copious pair produc- 
tion near the star. Since the plasma properties 
near the light cylinder seem to be insensitive to 
how exactly the near-star avalanche occurs, AE 
can be used directly to predict the lightcurves 
and spectra. It is already known (Bai, Spitkovsky 
2010) that postulating gamma-ray emission from 
the current layer beyond the light cylinder gives 
the lightcurves crudely consistent with observa- 
tions. 

4.2. Millisecond pulsars 

Estimating pair production near the light cylin- 
der for millisecond pulsars in the same manner as 
above, we also classify them as weak. But the 
near-star avalanche in millisecond pulsars prob- 
ably requires careful modeling, because by the 
usual criteria (Ruderman, Sutherland 1975) the 
avalanche is barely going. The near-star avalanche 
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can noticeably contribute to the observed emis- 
sion (although the current layer emission cannot 
be subdominant). 

4.3. Strong Pulsars 

Below we describe one scenario for what might 
be happening near the light cylinder in strong pul- 
sars - the "light-cylinder avalanche" . In this sce- 
nario, careful modeling of kinetics seems to be a 
must for predicting the emission. 

For the sake of the argument, take the light 
cylinder radius R ~ 3 x 10 8 cm and the mag- 
netic field near the light cylinder B ~ lOOkG. 
This would be a pulsar with L s d ~ 3 x 10 37 erg/s. 
Let e\L s d be the primary gamma-ray luminosity 
at about lGeV. Neglecting geometry, we estimate 
the lGeV gamma-ray density at the light cylinder 



47ri? 2 c(lGeV) 



10 12 eicm" 3 . 



(18) 
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Now send a test IkeV X-ray photon into the 
light cylinder region. With a probability 



P 



3 x 10" 5 ei 



(19) 



the test X-ray photon will produce a pair - an 
electron and a positron of ~GeV energy. In the 
lOOkG magnetic field, the pair will emit N ~ 10 6 
X-ray photons of ~keV energy . If pN > 1, we get 
an avalanche. 

In the simplest model, the pulsar balances at 
the threshold of the light-cylinder avalanche, e ~ 
e c ~ 0.03. For e > e c , copious pair production near 
the light cylinder opens the field lines (Fig. 2.) and 
the efficiency e drops. For e < e c , the avalanche 
disappears, pairs are no longer produced near the 
light cylinder and the field tries to return to the 
high-efficiency weak-pulsar configuration (Fig.l). 
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A. Numerical AE in axial symmetry 

The code used to generate the figures can be found at |http://cosmo. nyu.edu/an drei/AE/ 
We numerically integrate the AE equations ([TJ El [5J [HEH} with diffusive regularizations. It would be 
good to have a code which gets rid of the artificial diffusivities (3 and D (see below). 

A.l. Variables and Fields 

We assume axial symmetry. The coordinates are 

• t - time 

• (r, 9, z) - cylindrical coordinates 

• R = ^/r 2 + z 2 - spherical radius 

The fields arc 

• p± - positron/electron density (charge normalized), defined for R > r s 

• j± - positron/electron current, defined for R > r s 

• j - electric current 

• E, B - EM field 

A. 2. Parameters 

The parameters of the model arc 

• r s - the radius of the star 

• - the angular velocity of the star 

• a s - the electrical conductivity of the star 

• j e - permanent external current in the star (in the code, we use a purely toroidal current which scales 
linearly with r) 

• ai, oli — the pair production rates near the star and at the light cylinder 

• E c - the critical electric field for the inner pair production 

• 13 - "toroidal diffusivity" of the EM field 

• D - diffusivity of charges (in the code, we use diffusivity which scales linearly with R) 

A. 3. Equations 

The AE equations read 



G 



d t p± + r 1 d r (rj± r ) + d z {j± z ) = Q, R> r 



S, 



d t B r = d z Eg, 

d t B e = -d z E r + 8 r E z + (r- l d r {rd r Be) + d 2 z B e - r~ 2 Bg) , 
8 t B z = -r^drirEe), 

d t E r = -d z B g - j r , 

d t E 9 = d z B r - d r B z - jo + /3 {r- l d r {rd r Ee) + d 2 z Eg - r^Eg) , 
d t E z =r- 1 d r (rB e )-j z , 

j =j e +o s (E + VLr6 x B) , R<r s , 

j=j+-j-, _ R>r s , 
j± = P±v± - DVp±, 

v± = ExB±(Bg+^E) ; B 2_ E ^ = S 2_ E 2 t B E = B ■ V, E > 0, 

Q = ctiQi + ol^Qi-, 

Qi = 6{lAr s - R)9(E - E c )(l - ||), 
Q2 



(Al) 



En 

(SIR) 2 
l + (S1R) e ■ 



Note: at high a s , one can use a non-relativistic Ohm's law inside the star. For us the star is just a black 
box of high conductivity, and the detailed form of conductivity is irrelevant. 

A. 4. Initial and Boundary Conditions 

At t = all fields are equal to zero. 

We solve the problem inside a cylinder r < a, \z\ < b. The boundary conditions are absorbing for the 
density and outgoing for the EM field 

p± = 0, R = r s , 
p± = 0, r = a, 
P±=0, \z\=b, 

E r = -B e , z = -b, 
Eg = B r , z = -b, 

E r = B e , z = b, 
Eg = -B r , z = b, 

Eg = B z , r = a, 
E z = -Bg, r = a, 

Note: because of the diffusive term in the Bg equation, (Al), (A2) is not a correctly posed problem. 
But since (3 is just a regularization, we coded (Al), (A2) in a way which required no (explicit) additional 
boundary conditions. 



(A2) 
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